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Abstract 

An approach is presented for carrying out 
reliability-based design of a metallic, plate-like wing to 
meet strength and flutter requirements that are given in 
terms of risk/reliability. The design problem is to 
determine the thickness distribution such that wing 
weight is a minimum and the probability of failure is 
less than a specified value. Failure is assumed to occur 
if either the flutter speed is less than a specified allow- 
able or the stress caused by a pressure loading is greater 
than a specified allowable. Four uncertain quantities 
are considered: wing thickness, calculated flutter speed, 
allowable stress, and magnitude of a uniform pressure 
load. 

The reliability-based design optimization approach 
described herein starts with a design obtained using 
conventional deterministic design optimization with 
margins on the allowables. Reliability is calculated 
using Monte Carlo simulation with response surfaces 


^Aerospace Engineer, Analytical and Computational Methods 
Branch, Senior Member, AIAA 

^Aerospace Engineer, Analytical and Computational Methods 
Branch, Member AIAA 

^Aerospace Engineer, Aerodynamics, Structures, and Materials 
Department 

Copyright © 2001 by the American Institute of Aeronautics and As- 
tronautics, Inc. No copyright is asserted in the United States under 
Title 17, U.S. Code. The U.S. Government has a royalty-free license 
to exercise all rights under the copyright claimed herein for Govern- 
mental Purposes. All other rights are reserved by the copyright owner. 


that provide values of stresses and flutter speed. During 
the reliability-based design optimization, the response 
surfaces and move limits are coordinated to ensure 
accuracy of the response surfaces. 

Studies carried out in the paper show the relation- 
ship between reliability and weight and indicate that, 
for the design problem considered, increases in reliabil- 
ity can be obtained with modest increases in weight. 

Introduction 

Advanced aerospace vehicles of the future will be 
built to meet design requirements that are defined in 
terms of the likelihood, or probability, that the vehicle 
will perform as intended. In this paper, that likelihood 
is referred to as the reliability, and the associated design 
approach is referred to as reliability-based design. 

Reliability-based design has several advantages 
over traditional deterministic design. For example, 
with these methods, measures of reliability and safety 
are available during the design process and for the final 
design. This information allows the designer to pro- 
duce a consistent level of safety and efficiency 
throughout the vehicle — no unnecessary over-designs 
in some areas. As a result, designers may be able to 
save weight while maintaining adequate reliability and 
safety. In addition, in reliability-based design it is pos- 
sible to determine the sensitivity of the reliability to 
design changes that can be linked to changes in cost. 
Designers could carry out trade studies involving 
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reliability, weight, and cost. Also, for the same cost, it 
may be possible for aerospace vehicles to be made 
more reliable and safer than with traditional design 
approaches, or, for the same safety and reliability, the 
vehicle may be made at a lower cost. The primaiy dis- 
advantages of reliability-based design are that a sub- 
stantial amount of both data and computations are 
needed. 

To address these disadvantages, this paper provides 
an example of data collection and illustrates computa- 
tional techniques that help make reliability-based 
design tractable. In addition, this paper provides insight 
into reliability-based design for an example design 
problem involving two disciplines— aeroelasticity and 
stress analysis. 

This paper presents a design procedure by which a 
plate-like wing with variable thickness can be designed 
to be minimum weight and to meet both a subsonic 
flutter requirement and a separate strength requirement. 
However, unlike earlier, deterministic work (e.g., 
Refs. 1-4), this design procedure considers several 
uncertainties and accounts for them by using a prob- 
abilistic approach. Also, instead of an inequality con- 
straint on failure (flutter speed too low and/or stress too 
high), the inequality constraint is on the probability of 
failure. The design variables define the wing thickness 
distribution. 

First, the paper provides a description of the design 
problem. Then the details of the approaches used and 
the results of the studies are presented. The studies 
show the relationship between reliability and weight 
and indicate that, for the design problem considered, 
increases in reliability can be obtained with modest 
increases in weight. 

Description of Wing Design Problem 

The objective of the study is to obtain the 
minimum-weight plate-like wing that satisfies design 
requirements on both flutter and strength. This exam- 
ple was chosen because a similar metal plate can serve 
as the center portion, and primary stiffness, of a wind 
tunnel model. The appropriate aiifoil shape can be cre- 
ated by bonding balsa wood or a foam material to the 
metal plate and tooling the combination. The thickness 
distribution of the metal plate can be tailored so that the 
complete model has the desired stiffness distribution. 

The wing planform is shown in Figure 1. The wing 
semispan is 60 in., the wing root is 24 in., the tip is 
12 in., and the sweep of the leading edge is 14 degrees. 


The wing thickness distribution is taken to be quadratic 
in both the chordwise and spanwise directions. The 
quadratic thickness distribution is defined by the thick- 
nesses at the nine locations shown in Figure 2. Loca- 
tions 4, 5, and 6 are at mid semispan. Locations 2, 5, 
and 8 are at the mid points of the local chord. These 
nine thicknesses are the design variables. The expres- 
sion for the thickness distribution is presented and 
explained in Appendix A. The minimum gage require- 
ment for the thickness is 0.150 in. 



wing. 



Figure 2. Location of nine thicknesses that define the 
thickness distribution of the plate-like wing. 

There is one flutter speed requirement and one 
strength requirement. The flutter speed requirement is 
subsonic. The strength requirement is applied to each 
finite element. The loading for the strength analysis is 
uniform pressure over the surface of the wing. The 
wing planform is fixed; design variables control the 
thickness distribution. There are four uncertainties: 
wing thickness, calculated flutter speed, allowable 
stress, and magnitude of the uniform pressure load. 
Instead of an inequality constraint on failure (flutter 
speed too low and/or stress too high), the inequality 
constraint is on the probability of failure. 

The computer code MSC/NASTRAN was used to 
calculate both the stresses (Ref. 5) and the flutter speed 
(Ref. 6). The finite element grid was obtained by 
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dividing the wing into 18 equally-spaced segments 
spanwise and 9 segments choidwise — making a total of 
162 finite elements. The finite element grid is shown in 
Figure 3. The finite element that was used in the stress 
analysis is denoted CQUAD8. 

For the flutter analysis, the KE method, which is a 
modified K method, was used. Eight flutter modes 
were calculated. The frequency and damping for each 
mode were calculated at 99 reduced frequencies. The 
flutter speed for each of these modes was obtained from 
the reduced frequency at which the damping becomes 
negative. The Doublet-Lattice method was used to 
calculate the required aerodynamic matrices. The aero- 
dynamic grid was obtained by dividing the wing into 
10 equally-spaced segments spanwise and 5 segments 
chordwise. The aerodynamic grid is shown in Figure 4. 

Values of non-random quantities that were used in 
the wing design study are given in Table 1. These 
quantities were used in the stress analysis and in the 
flutter analysis. Definitions of the four types of random 


quantities that were used in canying out the probabilis- 
tic wing design study are given in Figure 5. The nine 
random thicknesses t;( random) were obtained by taking 
the input values of t, and adding a random increment 
from the uniform distribution shown in Figure 5. Note 
that the mean value of the increment is not zero. 
Further discussion of randomness in t; is given in 
Appendix B. 

The fourth type of random quantity defined in Fig- 
ure 5 is the calculated flutter speed, which depends on 
the thickness distribution (defined by t,( random)), on 
the aerodynamic grid, and on the flutter analysis 
method (e g., KE and PK). To account for the uncer- 
tainties associated with the aerodynamic grid and 
analysis method, an increment from a uniform distribu- 
tion was added to the flutter speed calculated using 
tj( random) The increment was based on the assump- 
tion that the flutter speed calculated with other aerody- 
namic grids and by other methods is within ±3% of the 
flutter speed calculated with the grid and method used 
here. 


Table 1. Values of the non-random quantities that define the wing design problem 


Quantities associated with wing metallic material 


Modulus of elasticity 

10 x 10*> psi 

Density 

0.10 lb/in 3 

Poisson’s ratio 

0.30 

Quantities associated with flutter analysis 


Reference air density 

0.0023 slug/ft3 

Mach Number 

0.45 

Reference Chord 

24 in. 

Structural damping (g) 

(over the range of frequencies from 0-1000 Hz) 

0.01 
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A 



Uniform Pressure Load 

Normal distribution 

Mean = 1 .0 psi 

Standard deviation = 0,05 psi 

Allowable Stress 

Normal distribution 

Mean = 50,000 psi 

Standard deviation = 2,500 psi 

Thicknesses 

/{random) = / ; ( given) + random increment from 
uniform distribution below 

A 


I 1 1 p. 

0.0025 +0.0055 

Calculated Flutter Speed 

Use /{random) to calculate flutter speed Vf. 
Fj(random) = Vf + random increment from 
uniform distribution below 

A 


1 1 ► 

— 0.03 V t + 0.03 V f 


Figure 5. Definitions of the four random quantities 
used in the probabilistic design problem. 

With the above information, the wing design 
problem can be stated as follows. Find the values of the 
thicknesses t; (i=l, . . . , 9) such that wing weight is a 
minimum and the probability of system failure (Pf) is 
less than a specified value (e g., 0.01). The objective 
function is the wing weight calculated using t, rather 
than the mean value of the wing weight calculated 
using tj( random). System fail uic is assumed to occur if 
any constraint is violated. The constraints are (1) the 
flutter speed must be greater than a specified allowable 


and (2) the von Mises stress at the centroid of any 
element must be less than the specified allowable. The 
design problem is an example of a series system in 
which there are many failure modes, any of which can 
cause the system to fail. 

Design Approach 

The design approach was a combination of 
an analyzer (MSC/NASTRAN), an optimizer (DOT, 
Ref. 7), and FORTRAN and C code for approximate 
analysis and probabilistic analysis. First, the method 
for selecting the starting point for the reliability-based 
optimization is discussed. Next, the analyses used to 
perform the reliability-based design are described. 
Then, possible sources of error and methods for error 
alleviation in this approach are discussed. 

The first step was to cany out a deterministic opti- 
mization with margins on the allowables. The allow- 
ables for both the deterministic and probabilistic 
optimization are shown in Table 2. Standard methods 
were used to cany out the deterministic optimization. 


Table 2. Allowables for both deterministic 
and probabilistic optimization 


Allowable 

Type of optimum 

Deterministic 

Probabilistic 

Stress 

44,000 psi 

50,000 psi 
(Mean) 

Flutter Speed 

550 ft/sec 

540 ft/sec 


The deterministic optimum served as a starting 
point for the probabilistic optimization. Monte-Carlo 
simulation (MCS) was used to calculate the reliability. 
Approximate analysis techniques in the form of 
response surfaces were used to provide estimates of the 
flutter speed and the stresses. The stresses were the 
von Mises stresses evaluated at the centroid of each of 
162 finite elements. A response surface was generated 
for this stress in each element. Therefore, with the 
response surface for the flutter speed, a total of 163 
response surfaces were used during the calculation of 
the reliability. (The type of response surface used is 
discussed in Appendix C.) In calculating the reliability, 
MCS received all input from the response surfaces 
rather than from MSC/NASTRAN. 

Two sources of error, MCS and the response 
surface calculations, affected the calculation of the 
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Figure 6. Move limit strategy during probabilistic 
optimization. 

reliability. In the following two paragraphs, these two 
sources of error are discussed along with the approach 
taken to minimize their effect 

In MCS, the reliability is computed from an integer 
count of the number of simulations that meet the reli- 
ability requirements. A large number of simulations are 
more likely to accurately reflect the distribution of the 
underlying random variables. As a result, a large num- 
ber of simulations ate more likely to provide an accu- 
rate prediction of the reliability. In addition, the 
optimizer requires derivatives of the reliability with 
respect to the design variables. These derivatives are 
computed by central differences using small changes in 
the design variables. Small changes in the design vari- 
ables may not change the count in an integer-based reli- 
ability calculation unless a large number of simulations 
are performed. In this paper, to maintain the required 
accuracy, one million MCS trials were used for each 
estimate of the reliability. 5 

In generating a response surface, points are 
selected to provide an accurate estimate of the response 
in the region of interest. Initially, several sets of 
response surfaces with different sampling region sizes 
were generated to study the accuracy of the response 
surfaces. Errors in the response surfaces were too large 
to use the entire design space for a sampling region. 
Because the error in the response surfaces must be 


ft is possible to estimate the error associated with a Monte Carlo 
simulation. For example, if a Monte Carlo simulation with one mil- 
lion trials indicates a reliability of 0.9900, there is a 95% probability 
that the true answer is between 0.9898 and 0.9902 (refs. 8,9). 


much less than the variation in the responses due to 
uncertainty, it was necessaty to use response surface 
sampling zones no larger than 0.01 inches for this 
problem. 

For this design problem, multiple response surfaces 
are used. When the optimizer selects a new design, the 
move limits are adjusted to match the region over 
which the response surfaces are accurate. The move 
limit strategy is illustrated in Figure 6, which shows a 
hypothetical optimization search in 2-D design space. 
Once the move limits are reached, new response sur- 
faces are generated, a new move direction is calculated, 
and new move limits are set. 

As mentioned previously, minimum gage thickness 
was chosen to be 0.150 in. Before calling MSC/ 
NASTRAN to generate the response surfaces, the wing 
thickness was calculated at the comets of each finite 
element. If the thickness was found to be less than 
minimum gage, it was set equal to 0.150 in. 

Results and Discussion 

First, the deterministic optimization is discussed. 
Then, the probabilistic optimization is discussed. 
Results from the deterministic optimization served 
as a starting point for the probabilistic optimization. 

Deterministic optimization 

In this study, two different deterministic optimum 
designs were found with approximately the same 
weight. For each of these two designs, the weight, 
flutter speed, largest von Mises stress (and element 
number), and the design variables are presented in 
Table 3. 

Contour plots showing the thickness distributions 
for the two designs are presented in Figure 7. Substan- 
tial parts of both wings are at minimum gage (0.150 in). 
Both designs are thick along the leading edge near mid- 
semispan. In addition, design 1 is thick at the root near 
mid-chord, while design 2 is thick at the root near the 
trailing edge. Plots showing the von Mises stress at the 
centroid of each element are presented in Figure 8 . The 
stress distributions are similar, in many respects, to the 
thickness distributions. In design optimization such as 
this, the procedure adjusts thicknesses to simultane- 
ously improve both the local strength and the load 
paths. The load paths for designs 1 and 2 are different. 

To explore why there appear to be two converged 
deterministic optimum designs, a technique was used 
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Table 3. Quantities associated with deterministic optima 


Quantity 

Optimum 1 

Optimum 2 

Weight (lb) 

30.820 

30.818 

Flutter speed (ft/sec) 

550.06 

549.61 

Largest von Mises stress (psi) 

43,999. 


at centroid of element number. 

5 


Design variables (in) 



ti 

0.15000 

0.15000 

t2 

0.54880 

0.43030 

l 3 

0.15000 

0.51194 

ft 

0.60145 

0.56690 

t5 

0.19433 

0.21366 

t6 

0.15000 

0.15000 

Q 

0.15000 

015000 

*8 

0.15000 

0.15000 

k) 

0.15000 

0.15000 



+ ^ < oo m c-1 O-v co © t> Tt -h oo m 

£•-. in «/"» ^ ^ cn m cn cn c-t <n cn 1 — ' 1 — 1 

^ 0000000 00 00000 

© 

Thickness, in. 

Figure 7. Thickness distaibution for optimum determi- 
nistic designs 1 and 2. 

that is similar to a technique used in Reference 10. A 
plane containing three designs— the two deterministic 
optimum designs plus a third design - is passed through 
design space. For the approach used herein, the third 
design has design variables that are two percent larger 



v-> rt o t— u-i <n o i ■ ir-, ci o r~ >/-i rt 

oi on r-i cK NC H On no On no 

-3- -3- co cn cn ol ot n) -h —h -h 

Stress, ksi 

Figure 8. Von Mises stress at centroid of each finite 
element for optimum deterministic designs 1 and 2. 

than the average of the corresponding design variables 
of the two deterministic optimum designs. On this 
plane, weight contours and constraint boundaries can be 
superimposed. Such a plot provides insight into various 
interactions in design space. 
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(a) Constraints are shown as a system constraint 
boundary that incorporates all constraints. 


(b) Several individual constraints are shown. Weight 
contours are not shown. 


Figure 9. Plane through design space with optimum deterministic designs 1 and 2. weight contours, and 
constraints. 


Such a chart is presented in Figure 9. In Figure 9a 
the weight is indicated by the color. Optimum determi- 
nistic design 1 is at the lower left, optimum determinis- 
tic design 2 is at the lower right, and the third design is 
at the top of the trapezoid. The curve denoted “System 
constraint boundary” is the envelope of the flutter con- 
straint and the stress constraint in all elements. The 
combination of the weight contours, the constraint 
boundaries, and the fact that a straight line drawn 
between the two designs passes through infeasible 
design space indicate that the two designs are distinct, 
converged optima. 

Figure 9b shows the individual constraints that 
form the envelope constraint boundary that is shown 
in Figure 9a. For simplicity, weight contours are not 
shown in Figure 9b. Optimum design 1 is bounded by 
rising weight, by the flutter constraint, and by the stress 
constraints from elements 5 and 46. Optimum design 2 


is bounded by rising weight, by the flutter constraint, 
and by stress constraints from elements 6, 27, and 46. 

Probabilistic optimization 

The reliabilities of the two optimum deterministic 
designs were calculated using the probabilistic allow- 
ables given in Table 2. For designs 1 and 2, the reli- 
abilities are 0.8967 and 0.8990, respectively. 

Probabilistic designs were obtained using, as start- 
ing points, both deterministic optima. The two starting 
points produced two different probabilistic designs; the 
probabilistic designs are similar to the original determi- 
nistic designs. The probabilistic designs are given the 
same designation as the deterministic design starting 
points — that is, optimum probabilistic design 1 was 
obtained using, as a starting point, optimum determi- 
nistic design 1. Optimum probabilistic design 1 
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Table 4. Quantities associated with optimum 
probabilistic design 1. System reliability 
requirement is 0.9990 


Quantity 

Value 

Nominal Weight, lb 

31.7659 

Mean Weight, lb 

31.8982 

Pf (stress) 

0.001002 

Pf (flutter) 

0.000006 

Pf (system) 

0.001008 

Design variables 

ty 

0.155908 

ty 

0.151994 

ty 

0.151399 

U 

0.150000 

ty 

0.579366 

ty 

0.150245 

ty 

0.150000 

ty 

0.594500 

ty 

0.207701 


was obtained for one reliability— 0.9990. Optimum 
probabilistic design 2 was obtained for four 
reliabilities— 0.9000, 0.9900, 0.9990, and 0.9999. For 
each design, the weight, the probability of various 
modes of failure, and the values of the design variables 
are given in Tables 4 and 5 for designs 1 and 2, respec- 
tively. Note that, in the tables, the quantity Pf (stress) is 


the probability of failure by exceeding the allowable 
stress in any element. (Optimum probabilistic designs 
could have been obtained for several values of the reli- 
ability for either probabilistic design 1 or probabilistic 
design 2. Design 2 was selected to demonstrate the 
technology in this paper.) 

In Tables 4 and 5 and elsewhere in this paper, the 
quantity “Nominal Weight” is the weight calculated 
using the deterministic thickness design variables t; 
rather than the random thicknesses t,( random). The 
quantity “Mean Weight” is the mean value of the 
weight calculated using the random thicknesses 
tj( random) Random thicknesses are always used in 
calculating the reliability. 

The nominal and mean weights of design 2 
(Table 5) are presented as a function of probability of 
system failure (Pf) and reliability in Figure 10. Five 
thousand Monte Carlo simulations were used to deter- 
mine the mean weight. Note that the horizontal scale is 
logarithmic with respect to Pf. Below the values of Pf 
are the corresponding values of the reliability. Plotted 
in this manner, the nominal and mean weights are al- 
most linear with respect to both Pf and to the reliability. 
The curves show the difference between the two 
weights— about 0.12 lb or about 0.38 percent. The 
curves also show that large increases in reliability can 
be obtained with relatively small increases in weight 


Table 5. Quantities associated with optimum probabilistic design 2 
for foui' system reliability requirements 


Quantity 

System Reliability Requirement 

0.9000 

0.9900 

0.9990 

0.9999 

Nominal Weight, lb 

30.7185 

31.2389 

31.7918 

32.3702 

Mean Weight, lb 

30.8500 

31.3715 

31.9277 

32.5073 

Pf (stress) 

0.051390 

0.009774 

0.000012 

0 000100 

Pf (flutter) 

0.051910 

0.000230 

0.000986 

0000000 

Pf (system) 

0.100001 

0.010003 

0.000998 

0.000100 

Design variables 
tt 

0.150000 

0.150054 

0.150592 

0.153004 

ty 

0.429735 

0.440164 

0.457556 

0.467194 

ty 

0.515641 

0.525466 

0.536452 

0.547373 

U 

0.572689 

0.568498 

0.555382 

0.563198 

ty 

0.207682 

0.217424 

0.226726 

0.232476 

ty 

0.150139 

0.150123 

0.154583 

0.157350 

ty 

0.167267 

0.157169 

0.150243 

0.150000 

ty 

0.150091 

0.150000 

0.150487 

0.150072 

ty 

0.152205 

0.150905 

0.151110 

0.150675 
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System design requirements 


Figure 10. Nominal and mean weights of optimum 
probabilistic design 2 as a function of system require- 
ment on probability of failure (Pf ) and reliability (R). 
R= 1 - P f . 
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Weight, lb 

Figure 11. Distribution of weights produced by 
tj(random) for design 2, reliability = 0.9990. 



Figure 12. Thickness distribution for optimum 
probabilistic designs 1 and 2. System reliability 
requirement is 0.9990. 


Contour plots that show the thickness distributions 
of probabilistic designs 1 and 2 for a reliability of 
0.9990 are presented in Figure 12. These thickness 
distributions were calculated using the deterministic 
design variables tj rather than the random thicknesses 
tj(random). If the same terminology that was used 
above to describe the weight is used here, the thickness 
distributions shown in Figure 12 could be referred to as 
Nominal Thicknesses rather than Mean Thicknesses. 
These thickness distributions are similar to those in 
Figure 7 for the deterministic designs. 


In Tables 4 and 5, six significant digits are pre- 
sented because that precision is needed to obtain a 
minimum weight design for a reliability of 0.9999. To 
be consistent, six significant digits are presented for all 
values of design reliability. 

The range of weights produced by tj(random) is 
indicated in Figure 11. These data are for design 2, 
system reliability = 0.9990. Five thousand Monte Cario 
simulations were used to determine the distribution. 
The nominal weight and the mean weight are also 
shown. The range of possible weights is greater than 
the difference between the nominal and mean weights. 


Figure 13 provides detailed information about the 
stresses in probabilistic designs 1 and 2 that were 
designed for a reliability of 0.9990. This figure indi- 
cates the probability that the von Mises stress exceeds 
the allowable stress at the centroid of each finite 
element. Whereas Tables 4 and 5 give the probability 
of exceeding the allowable stress in any element. 
Figure 13 gives the probability of failure for each ele- 
ment. The results show that there are only a few ele- 
ments where the stress has a relatively high likelihood 
of exceeding the allowable stress. These elements are 
the same as those that provide the active constraints for 
the deterministic designs. 
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Figure 13. Probability that von Mises stress exceeds 
the allowable stress at the centroid of each finite ele- 
ment. Designs are optimum probabilistic designs 1 and 
2 for system reliability requirement of 0.9990. 
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Figure 14. Plane through design space with optimum 
probabilistic designs 1 and 2, weight contours, and 
curves of constant system reliability. 


Finally, a chart similar to Figure 9 is shown in Fig- 
ure 14 to help explain the differences between prob- 
abilistic designs 1 and 2 for a reliability of 0.9990. The 
third design in Figure 14 was obtained the same way it 
was obtained for Figure 9. The difference between 
Figure 9 and Figure 14 is that in Figure 14 the con- 
straint curves are system reliability. Although the 
design space in Figure 14 is not exactly the same as in 
Figure 9, it is close enough to see how various failure 
modes affect the reliability. For example, the sharp 
changes in the reliability in the right portion of Fig- 
ure 14 are caused by the multiple stress failure modes 
shown near design 2 in Figure 9b. 

Concluding Remarks 

An approach for carrying out reliability-based 
design of a metal, plate-like wing that can serve as the 
center portion and primary stiffness of a wind tunnel 
model is presented. The plate is designed to meet 
strength and flutter requirements defined in terms of 
risk/reliability. Although this study does not address all 
the issues associated with designing this type of wind 
tunnel model, this study can serve as the first step in 
developing a capability to account for uncertainties in 
the design of such a wind tunnel model. The design 
approach can also be used for preliminary design of 
wing structures for aircraft with thin wings. 

The design problem that was studied can be 
defined as follows. Determine the thickness distribu- 
tion such that wing weight is a minimum and the prob- 
ability of failure is less than a specified value. Failure 
is assumed to occur if either the flutter speed is less 
than a specified allowable or the stress caused by a 
pressure loading is greater than a specified allowable. 
There are four uncertainties: wing thickness, calculated 
flutter speed, allowable stress, and magnitude of a uni- 
form pressure load. Data for uncertainties in the thick- 
ness distribution of the metal wing were obtained from 
a NASA machine shop that has made such models. 

The key steps of the design approach are as fol- 
lows. The reliability-based design optimization started 
with a design obtained using conventional deterministic 
design optimization with margins on the allowables. 
The reliability was calculated using Monte Cario simu- 
lation (MCS) with one million trials. MCS interacted 
with response surfaces that provided values of stresses 
and flutter speed. During the reliability-based design 
optimization, the response surfaces and move limits 
were coordinated to ensure accuracy of the response 
surfaces. 
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The studies produced two reliability-based designs 
that met the design requirements and had nearly the 
same weight. A novel graphical approach is used to 
show how the two designs are distinct, converged 
optima. 

The studies showed the relationship between reli- 
ability and weight and showed that, for the design 
problem considered, large increases in reliability could 
be obtained with relatively small increases in weight. 

Appendix A. Mathematical Expression that Represents 
Wing Thickness Distribution 

The thickness varies over the wing as a biquadratic 
polynomial in x and y as follows 

2 2 

t = a\ + ^ 2 ^ + ^ 3 xy + a^y + a^x +a^x y 
2 2 2 2 

+ dqxy +a$y +a<^x y (Al) 

The unknown coefficients cij (i = 1,9) in Equa- 
tion (Al) can be determined from thickness values at 
nine locations in the wing. For the present study, the 
nine locations are shown in Figure 2 and in Fig- 
ure Al(a). Using Lagrange finite element shape func- 
tions and the thickness at these nine locations on the 
wing, the thickness distribution r(§,q) can be described 
in a curvilinear coordinate system (Fig. Al(b)) as 

t (I- rj) = (c x + c 2 q + c 3 rf j + (c 4 + c 5 q + c 6 i] 2 )§ 

+ (c 7 +c 8 rj + c 9 rj 2 )§ 2 (A2) 



(a) Wing planform in Cartesian (x, y) 
coordinate system. 


where 


C 1 = f 5 

(A3) 

c 2 = -{h- t i) 

(A4) 


(A5) 

c 4 = —{ t 6~ t 4) 

(A6) 

c 5 = —(ti + t 9 -t 3 -t 7 ) 

(A7) 

c 6 = —{U ~ %) + —{h + f 9 _ h ~ f l) 

(A8) 

C 1 =^( ? 4 +t 6)~ t 5 

(A9) 

c S = -^{ t 2~ t 8) + '^{ t l + t 9~ l l~ h) 

(A10) 

c 9 = { 5 "( ? 2 +t 4 +t 6 + tg) 
+ (r( r i +t 3 + t 7 + t 9 ) 

(All) 


— 1—1 — 1,1 



(b) Location of nodes in curvilinear 

(|, rf) system. 


Figure Al. Wing planform and Lagrange isoparametric transformation. 
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The equation relating (x,y) and (f,j ;) can be written 
as 


V = 



1 = 


2L(x - y tanfl) 
c r L-(c r -c t )y 


1 


(A 12) 

(A 13) 


The quantities L, 0, c r , and ct in Equations (A 12) and 
(A13) are shown in Figure Al(a). 

When response surfaces were created for stress and 
flutter, the wing was modeled using 162, eight-node 
quadrilateral finite elements. The finite element analy- 
sis code MSC/NASTRAN requires the thickness at 
each of the four comer nodes of each quadrilateral 
finite element as input. The thicknesses of the comer 
nodes were obtained using Equation (A4) for a given 
set of thickness (fj to t 9 ) . However, if a calculated 
thickness is less than the minimum thickness 0.15 in., 
the thickness at that node is set equal to 0.15 in. 

Appendix B. Origin of Randomness in 
Thickness— Machine Shop Measurements 


Figure B1 presents a histogram of the differences 
between a measured dimension and the specified 
dimension. Positive numbers indicate an excess of 
material. The results show that the number of points 
with excess material, and the magnitude of the 
excesses, exceeds the number and magnitude of the 
undercuts. This distribution reflects the approach used 
to produce the model. After the model was machined 
and carefully measured, the model and measurements 
were shown to the customer. The customer could 
accept the job, or could request some rework. Since it 
was much easier to rework the job by removing mate- 
rial than by adding material, it was much more likely 
that the machinist would leave excess material than 
undercut the job. Therefore, the distribution shown in 
Figure B1 is not symmetric about a difference of zero. 

Regarding the undercut data at the left in Fig- 
ure Bl, since the customer specified a tolerance of 
±0.0025 in., only a small number of occurrences fell 
outside that tolerance. The data in Figure Bl was 
approximated by a uniform distribution between 
-0.0025 and +0.0055, which is shown in Figure 5. 

Appendix C. Response Surface Used in Calculations 


Data describing the uncertainties associated with 
machining a wind tunnel model were obtained from 
quality control measurements of a wind tunnel model 
that had been machined and measured at NASA 
Langley Research Center. These measurements were 
taken at 14,805 points and were compared with the 
plans provided by the customer, the test engineer. The 
plans called for a tolerance of ±0.0025 in. 



Difference between dimensions, in. 


In this paper, Monte Carlo simulation (MCS) is 
used to calculate the reliability. Unfortunately, MCS 
requires a large number of function evaluations to 
obtain the required accuracy. Because the stress and 
flutter analyses used in this paper arc computationally 
expensive, it would not be practical to perform exact 
analyses for eveiy simulation. Instead, an alternate 
method was used to calculate the stresses and flutter 
speed. That method was based on response surface 
approximations. 

In the response surface approach, approximate 
polynomial functions are fit to the response over the 
design space. MCS interacts with these relatively sim- 
ple functions rather than with the more complex exact 
analyses. In the present work, one quadratic polyno- 
mial is used to represent the stress response in each 
of the 162 finite elements, and one quadratic polyno- 
mial is used to represent the flutter speed, for a total of 
163 response surfaces. A quadratic polynomial with 
9 design variables (v) requires k = (v + l)(v + 2)12 = 55 
coefficients. 


Figure Bl. Distribution of difference between 
measured dimensions and a specified dimension 
for a wind tunnel model. Difference is given by 
measured dimension minus specified dimension. 


In this study, the 55 coefficients of the polynomial 
were estimated by using a face-centered central com- 
posite design approach (Ref. 11). In this approach, 
design points are established at the center of the design 


12 

American Institute of Aeronautics and Astronautics 



space, at the comers of the design space (the intersec- 
tion between each of the design variable limits), and at 
the face-centers of the design space (points determined 
by setting one variable to its upper or lower limit with 
all other variables set to the average of their upper and 
lower limits). The number of design points, p, required 
by this procedure is p = 2 V + 2v + 1 = 531 points. The 
method of least squares is used to solve for the 55 coef- 
ficients by minimizing the error between the 531 design 
points and the response polynomial. 
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